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Effect of uniaxial strain on the Drude weight of graphene 

F. M. D. Pellcgrino a < b ' c < d , G. G. N. Angilella a > b ' c < d > c * and R. Pucci a ' c 

a Dipartimento di Fisica e Astronomia, Universita di Catania, Via S. Sofia, 64, 1-95123 
Catania, Italy; h Scuola Superiore di Catania, Via S. Nullo, 5/i, 1-95123 Catania, Italy; 
C CNISM, UdR Catania, 1-95123 Catania, Italy; d CNR-IMM, Z.I. VIII Strada 5, 1-95121 
Catania, Italy; °INFN, Scz. Catania, 1-95123 Catania, Italy; 

{September 20, 2010) 

We study the dependence on the strength and orientation of applied uniaxial strain of 
the Drude weight in the conductivity of graphene. We find a nonmonotonic dependence on 
strain, which may be related to the proximity to several strain-induced electronic topological 
transitions. Given the relation between the Drude weight and the long-wavelength plasmon 
frequency, such a strain dependence can be evidenced by infrared spectroscopy measurements. 
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Introduction. The recent fabrication in the laboratory [l[ of graphene — an atom- 
ically thick single layer of carbon atoms in the sp 2 hybridization status — has 
spurred a wave of research activity both on the experimental and on the theoret- 
ical side. Graphene is characterized by a honeycomb lattice, made of two inter- 
penetrating triangular sublattices. This gives rise to unusual electronic properties. 
The band structure of graphene thus consists of two bands, touching at the Fermi 
level in a linear, cone-like fashion at the so-called Dirac points ±K. Indeed, the 
low-energy quasiparticles in graphene follow much the same physics as relativis- 
tic massless particles, thus allowing the observation of several quantum-relativistic 
electrodynamic effects in a solid state system 

Equally remarkable are graphene's outstanding elastic properties. Despite its 
quasi-two-dimensional character — that should in principle even prevent its very 
existence, owing to Mermin- Wagner's theorem 0] — graphene is characterized by 
exceptional tensile strength and stiffness 0], with a record breaking strength of 
~ 40 N/m 14], which nearly equals the theoretical limit. The Young modulus of 
~ 1.0 TPa [J] places graphene amongst the hardest materials known. Recent ah 
initio calculations as well as experiments [1] have demonstrated that graphene 
single layers can reversibly sustain elastic deformations as large as 20%. It has been 
even proposed that graphene-based nanodevices could be engineered, by tailoring 
the desired electronic or transport properties on the strain-induced modifications 
of the deformed graphene sheet (origami electronics) fiol ]. It has been predicted 
that uniaxial strain should influence the optical conductivity of graphene [H[ as 
well as its reflectivity in the optical range 121 ]. 

Here, we will be concerned on the strain effects on the Drude weight of the 
conductivity of graphene. This quantity can be connected by means of an effective 
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/-sum rule to the dis per sion relation of plasmons [lit ] , which has been studied also 
under applied strain [14]. 

Model. The band structure of electrons residing on a honeycomb lattice can be 
described in terms of the tight-binding Hamiltonian 



H = J]^a t (R)6(R + ^) + H.c., (1) 



where a'(R) is a creation operator on the position R of the A sublattice, b(R + 8e) 
is a destruction operator on a nearest neighbor (NN) site R + 8e, belonging to 
the B sublattice, and 81 are the vectors connecting a given site to its nearest 
neighbors, their relaxed (unstrained) components being 8^ = a(l,\/3)/2, 8 



2 



a(l, — \/3)/2, 8 3 = a(— 1,0), with a = 1.42 A, the equilibrium C-C distance in 
a graphene sheet 0. In Eq. ©, t e = t(S e ), £ = 1,2,3, is the hopping parameter 
between two NN sites. In the absence of strain they reduce to a single constant, 
U = t , with t = -2.8 eV 0. 
In terms of the strain tensor e = e(e,6) [l6[ the deformed lattice distances are 

related to the unstrained ones by 8^ = (D + e) • 8^ . Here, 9 denotes the angle along 
which the strain is applied, with respect to the x lattice coordinate, e is the strain 
modulus. The special values 9 = and 9 = ir/6 refer to strain along the armchair 
and zig zag directions, respectively. 

The resulting tight-binding band structure features two bands EkA (A = 1,2), 
which close to the Fermi level can be approximated by Dirac cones, touching at 



±ko [111 ]. While at e = these coincide with the high symmetry points ±K at 
the vertices of the hexagonal first Brillouin zone (1BZ) of graphene, their posi- 
tions move away from ±K with increasing strain. Depending on e and 9, they may 
eventually merge into the midpoints of one of the sides of the 1BZ. In that case, 
the cone approximation breaks down, and a finite gap opens at the Fermi level 
[III ]. F° r intermediate strains, the constant energy contours of E^x can be grouped 
according to their topology, and are divided by three separatrix lines [111 ] . This 
corresponds to having three distinct electronic topological transitions (ETT) [l7r - 
Hil ]. which are here tuned by uniaxial strain, as has been suggested to be the case 
in other quasi-2D materials, such as the cuprates [13, HH and the Bechgaard salts 
[12]. Correspondingly, while the density of states (DOS) is characterized by a linear 
behaviour close to the Fermi level, because of the linearity in at k = ±ko, a 
detailed analysis beyond the cone approximation reveals the presence of logarith- 
mic singularities in the DOS at higher energies, that may be connected with the 
proximity to the ETTs [HI ]. This behaviour is reproduced in the frequency depen- 
dence of the longitudinal optical conductivity a(ui) and of the optical reflectivity 
R(cj), which have been studied at increasing strain modulus and different strain 



directions [ll|, [12| . 



Results and conclusions. We are interested in the longitudinal response to an 
electric field forming an angle ip with the x crystallographic direction. This is given 
by the dynamical conductivity [H| 



• 9 -2 

ie n ic 



o- w (uj) = + — lim 11$ ' (q, -q, u), (2) 

moo uj <}->0 w 



where II^,(q, — q, u) is the retarded current-current (jj) correlation function, e 
and m are the electron charge and mass, respectively, and n is the carrier density. 
This is related to the density-density (pp) correlation function U pp through the 
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continuity equation. One obtains 

ie 2 oj 2 

a w (co) = — lim ^n^(q, -q, u). (3) 

Letting oj — > u + i0 + , and extracting the real part, one recognizes the Drude weight 
as 

to 2 

D w = Tie 2 lim -^-RelF'tq, -q,0). (4) 

<j->0 q A 

Within linear response theory, collective excitations are described by the analytical 
properties of the appropriate correlation function. At lowest order, this can be 
expressed as a polarization bubble. The random phase approximation (RPA) is 
then the simplest diagrammatic procedure to include infinite orders in the dielectric 
screening [23J. At low frequency and small wavevectors, the k- space integration 
involved in the definition of W p can be performed analytically [141 ] . In that limit, 
the Fermi line in general reduces to an ellipse of constant energy, /i. One finds 



D w = 4pa 



l /cos (<p — rj) sin (tp — rj) 



^ A c Pi To + 



A 2 B 2 



(5) 



where o~q = ne 2 /2h is the so-called universal interband electrical conductivity of 
neutral graphene [13], A c = 3\/3a 2 /2 is the area of the unit cell, p\ is the strain- 
dependent prefactor in the linear dependence of the DOS on the chemical potential 
at low energy, p(ji) = g s p\\ji\ [ill ], with g s = 2 for the spin degeneracy, and A and B 
denote the semiaxes of the constant energy ellipse, with A~ 2 = =(7 — \/ a 2 + (3 2 ), 
B- 2 = I( 7 + ^TW); a = -W(tj + t 2 - 2i§), (3 = ^a 2 (t 2 - t 2 ), 7 = 



2 ' 

n 

2\l ' V" 1 y h " — 2" 7"! 1 "2 ""'ill r — 2" 

\a 2 (t 2 + 4 + t 2 3 ); and cos(2r/) = \a\/y/a 2 + /3 2 , sin(2?7) = \a\fi / a 2 + (3 2 . 



2 l 

In Eq. ([5]), the quantity between square brackets goes to unity in the limit e — > 0, 
where in particular p\ = 2/(7r\/3t 2 ). From Eq. ([5]), it follows that D w attains its 
maximum values whenever tp — rj = n/2 (modulo 7r), and its minimum values 
whenever <p> — r\ = (modulo it). It turns out that rj = 9 (the angle defining 
the direction of applied strain) in the zig-zag and armchair cases, whereas rj ~ 9 
in the generic case. Whilst Eq. ([5]) has been derived within the (deformed) cone 
approximation^which implies low strain, experimental results as well as ab initio 
calculations [5|-]8|] have shown that both lattice and electronic structures of graphene 
are stable with respect to strains as high as 20 %. We therefore expect Eq. ([5]) to 
be robust with respect to applied strain, within such strain ranges. In particular, 
the main strain dependence is summarized by p\, which modulates the overall 
angular dependence of D w . Such a quantity has been shown to increase steadily, 



and eventually diverge close to the cone breakdown [ll( . Also the ellipse semiaxes A 
and B do depend on strain, whose role is that of increasing the ellipse anisotropy, 
with A = B when e = 0. Such a dependence of the Drude weight on applied 
uniaxial strain is amenable to experimental verification. In fact, the prefactor £b\ 
in the long- wavelength dispersion relation of low-ener gy p lasmons in graphene [l3 |. 
w q = Coi^/qa, is related to the Drude weight through [131 ] 

D w _ 27re e r 



where eo is the dielectric constant in vacuum, e r = (ei + £2)/2, with e±, €2 the 
relative dielectric constants of the two media sandwiching the graphene sheet. 
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In conclusion, we have found a nonmonotonic strain dependence of the Drude 
weight, related to the longitudinal polarization of graphene. This can be described 
as an effect of the proximity to several strain-induced ETTs. Given the connection 
between the Drude weight and the long- wavelength plasmon frequency 13l . Il4j ]. 
such a dependence could be evidenced by infrared spectroscopy [251 ] . 
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